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RESUMEN 
ABSTRACT 

, We present the multi-color, five-year light curves and the first radial velocities 

C/5 ■ of the near-contact binary system KR Cyg. We derived the masses of the 

components as 2.88±0.20 M and 1.26±0.07 M and the radii as 2.59±0.06 
Rq and 1.80±0.04 R©. Analyses of the UBVR light curves and the radial ve- 
locities indicate that none of the components exactly fill their corresponding 
Roche lobes. We have calculated the distance to the system of KRCyg as 
^ 411±12 pc using the observed apparent UBV magnitudes and the bolometric 

O f corrections for the component stars. We also searched for the empirical deter- 

mination of albedo and effective temperature of the cooler, less massive star 
of KR Cyg, and of two similar near contact binaries AK CMi, and DO Cas. 
The residuals between the observed and computed fluxes are attributed to 
the effect of mutual illumination which heats the surface layers of 
the illuminated star and does vary not only its bolometric albedo 
but also its limb-darkening coefficient and gravity-brightening ex- 
ponent. The analysis of the light curves shows that the effective albedos 
are generally smaller than that expected from an envelope of convective star, 
being mostly departed from the theoretical value at the B passband. As 
the reflected light diminishes the effective temperature and, therefore, the 
luminosity of the irradiated star increase. The observed bluer U-B colors dur- 
ing primary minimum are attributed to the effects of mutual irradiation and 
multiple scattering processes which may alter several characteristics of these 
systems. 

Key Words: stars: binaries: close — binaries: eclipsing — binaries: general 
— binaries: spectroscopic — stars: individual: KR Cyg, AK CMi, DO Cas 



1. INTRODUCTION 



The light variabi lity of KRCyg (HD 333645, V 



tDiiity oi rvKUyg (tiu d3d04t>, v = 9 m .19) was discovered 
by iSchnellerl (|1931allbT ldY The light c urves obtai ned by photograph ic obser- 
vations were published by iNekrasoval (|1945f ) and IWachmannl (|1948l ) and the 
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visual observations by Lassovzskv ( 1936h . Gaposchkin (1953) and Tsesevich 
(|1954f ). Both visual and photographic observations showed that the light 
variations were originated from the eclipses. The light curves were very sim- 
ilar to th ose of Algol-typ e binaries. The first photoelectric observations were 
made by IVetesnikl (119651 ). He obtained blue and yellow light curves of the 
system. The passbands used in the observations were very similar to those 
Johnson's B and V filters. A preliminary solution of the light curves leads him 
to the conclusion that the system consists of a B9 and an F5 main-sequence 
stars. He classified the system as a j3 Lyrae type by visual inspection of 
the light curves. He also called attention some peculiarities in the course 
of the light curve. Later o n Vetesnik's two-color light curves were analysed 
bv IWilson fc RafertJ (Il980h using a conte mporary method, i.e. the Wilson- 



Devinney code (j Wilson k, Devinnevlll971l ). henceforth WD. The main differ- 



ence between the two solutions is the primary star's radius which is fifteen 
per ce nt larger in the later th an that found in the former analysis. In ad- 
dition, IWilson fe RafertJ |l980l) use an effective temperature for the primary 
star corresponding to a B7 star. The observ ations in the primary m inimum of 
Vetesnik's light curves were re-analysed by lA-Naimiv et al.l (|1985l ) using the 
Fourier techniques in the frequency domain. They find that, in contrary to 
the older estimates, more massive star has smaller radius than its less massive 
companion. 



Sipahi 



Gulmenl (|2004l ) published three-color photometric observations 
and the resultant light curves of KRCyg. The light curves were analysed 
individually and simultaneously by the WD code. They determined a pho- 
tometric mass-ratio to be 0.43 and suggested that the cooler, less massive 
component is nearly filling its corresponding Roche lobe. Recently, one of us 
(jSipahil 120121 ) collected all available times of light minimum and studied or- 
bital period of the system. She arrived at a result that the differences between 
the times of observed and calculated (O-C) show a periodic oscillation with 
an amplitude of 0.001 days and a period of about 76 years. This periodic 
change was attributed to a hypothetical third component. 



Shawl (|1990l ) defined a group of close binary stars named near-contact 



binaries (NCB) in which both components fill or nearly fill their critical Roche 
lobes. KR Cyg is included in the list of NCBs. According to his definition 
the stars of the near-contact binaries are near enough to each other to have 
strong proximity effects like W UMa type systems but are not in contact. 
As it is known the mass and energy exchange are being occurred between 
the components in the contact systems. This energy exchange is resulted 
in both components to come to nearly the same effective temperature. The 
main difference between the near-contact and contact binaries is the non- 
existence of the energy exchange through the neck of material connecting 
the components in the near-contact systems. The non-existence of energy 
exchange in the case of NCB allows a large temperature difference between 
the components. Some NCBs show photometric asymmetries such that the 
system is brighter at the first quarter than that at the second one. This 
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TABLE 1 

THE COORDINATES, APPARENT VISUAL MAGNITUDES AND THE 
SPECTRAL TYPES OF THE STARS OBSERVED. 



Star 


Alpha 


Delta 


V (mag) 


Sp 


KR Cyg 


20 09 06 


30 33 01 


9.23 


AOV 


HD 191398 


20 08 39 


30 20 15 


9.01 


A0V 


HD 333664 


20 09 18 


30 13 39 


9.58 


AO 



asymmetry is attributed to the primary star being one hemisphere, facing 
at phase about 0.25, is hotter than that at phase 0.75, or some material is 
obscuring the surface of the primary visible at phase 0.75. Therefore the NCBs 
are assumed to be lying in key evolutionary states of close binary systems. 



Shawl (| 1994T ) renewed the list of NCBs in which over 130 systems are given. 



He proposes that the NCBs are precursors of the A-type W UMa contact 
binaries. 

The main aim of this study is to analyse the 1999 and 2000 BVR, and 2003, 
2004, and 2005 UBVR light curves and the first radial velocities obtained by 
us. We derive the absolute physical parameters of the components and discuss 
light curve peculiarities by taking into account the mutual irradiation and 
therefore the changing in the value of albedo of the irradiated star for four 
passbands. 



2. OBSERVATIONS 

2.1. Photometric observations 

The observations were made with the High-Speed Three-Channel Pho- 
tometer attached to the 48 cm Cassegrain telescope at Ege University Obser- 
vatory HD 191398 and HD 333664, respectively, are taken as the comparison 
and check stars. These stars are nearly the same apparent magnitude and 
spectral type with the variable star. Moreover, they are very closely located 
with the variable on the plane of sky which make the effect of the atmo- 
spheric extinction on the differential magnitudes almost negligible. In Table 1 
the coordinates, apparent visual magnitudes and spectral types of the variable 
and the comparison stars are given. The observations were made using the 
BVR passbands in 1999 and 2000, and UBVR passbands in 2003, 2004 and 
2005. Though the comparison star is so close to the variable, the differential 
magnitudes, in the sense variable minus comparison, are corrected for atmo- 
spheric extinction. The nightly extinction coefficients were derived from the 
observations of the comparison stars in each passband. 



2.2. Spectroscopic observations 

Optical spectroscopic observations of KR Cyg were obtained with the Turk- 
ish Faint Object Spectrograph Camera (TFOSC) attached to the 1.5 m tele- 
scope on 4 nights (July- August, 2011) under good seeing conditions. Further 
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details on the telescope and the spectrograph can be found at http://www.tug.tubitak.gov.tr 
The wavelength coverage of each spectrum was 4000-9000 A in 12 orders, with 
a resolving power of A/ A A 7 000 at 6563 A and an average signal-to-noise ratio 
(S/N) was ~120. We also obtained a high S/N spectrum of the a Lyr (AO 
V) and HP 27962 fo r use as templates in derivation of the radial velocities 



( Nidever et alj|2002 ) 



The electronic bias was removed from each image and we used the 'crrc- 
ject' option for cosmic ray removal. Thus, the resulting spectra were largely 
cleaned from the cosmic rays. The echelle spectra were extracted and wave- 
length calibrated by u sing Fe-Ar lamp source with help of the IRAF echelle 
|Tonrv fc Davislll979l) package. 



The stability of the instrument was checked by cross correlating the spectra 
of the standard star against each other using the FXCOR task in IRAF. The 
standard deviation of the d ifferences between th e velocities measured using 



FXCOR and the velocities in lNidever et al.l (|2002l) was about 1.1 kms 1 



3. ANALYSES 
3.1. Radial velocities 

To derive the radial velocities for the components of binary system, the 
four TFOSC spectra of the eclipsing binary, obtained near the quadratures, 
were cross-correlated against the spectrum of a Lyr on an order-by-order basis 
using the FXCOR package in IRAF. The majority of the spectra showed two 
distinct cross-correlation peaks in the quadrature, one for each component of 
the binary. Thus, both peaks were fitted independently in the quadrature 
with a Gaussian profile to measure the velocity and errors of the individual 
components. If the two peaks appear blended, a double Gaussian was applied 
to the combined profile using de-blend function in the task. For each of the four 
observations we, then, determined a weighted-average radial velocity for each 
star from all orders without significant contamination by telluric absorption 
features. Here we used as weights the inverse of the variance of the radial 
velocity measurements in each order, as reported by FXCOR. 

We adopted a two — Gaussian fit algorithm to resolve cross-correlation 
peaks near the first and second quadratures when spectral lines are visible 
separately. Fig. 1 shows examples of cros-correlations obtained by using the 
largest FWHM at nearly first and second quadratures. The two peaks, cor- 
respond to each component of KRCyg. The stronger peaks in each CCF 
correspond to the more luminous component which has a larger weight into 
the observed spectrum. 

The heliocentric RVs for the primary (V p ) and the secondary (V s ) com- 
ponents are listed in Table 2, along with the dates of observation and t he cor 



respon ding orbital phases computed with the new ephemeris given by ISipahi 



(|2012l ). The velocities in this table have been corrected to the heliocentric 
reference system by adopting a radial velocity of 9.5 kms -1 for the template 
star a Lyr. The RVs listed in Table 2 are the weighted averages of the values 
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Fig. 1. Sample of CCFs between KR Cyg and the RV template spectrum around 
the first and second quadrature. 



TABLE 2 

HELIOCENTRIC RADIAL VELOCITIES OF 
KR CYGNI. THE COLUMNS GIVE THE 
HELIOCENTRIC JULIAN DATE, THE ORBITAL 
PHASE, THE RADIAL VELOCITIES OF THE TWO 
COMPONENTS WITH THE CORRESPONDING 
STANDARD DEVIATIONS. 

HJD 2400000+ Phase Star 1 Star 2 







v P 




V s 


a 


55751.36069 


0.8080 


73.1 


11.8 


-290.8 


11.3 


55796.38592 


0.0826 


-92.3 


7.8 


77.7 


11.2 


55835.36321 


0.2013 


-150.4 


6.8 


188.3 


9.5 


55835.46128 


0.3173 


-144.6 


9.9 


184.1 


9.9 



G 
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Fig. 2. Radial velocities of the components. 



obtained from the cross-correlation of orders #4#, #5#, #6# and #7# of 
the target spectra with the corresponding order of the standard star spec- 
trum. The weight Wi = 1/of has been given to each measurement. The 
standard errors of the weighted means have been calculated on the basis of 
the errors (o~j ) in the RV values for each order according to the usual formula 



(e.g lTopping I (|1972h ). The <n values are computed by fxc or according to 



the fitted peak height, as described by Tonrv fc Davis ( 19791 ) 



The measured radial velocities of both stars are plotted in Fig. 2 against 
the orbital phase. First we analysed the radial velocities for the initial or- 
bital parameters. We used the orbital period held fixed and computed the 
eccentricity of the orbit, systemic velocity and semi- amplitudes of the RVs. 
The results of the analysis are as follows: e=0.001±0.001, i.e. formally 
consistent with a circular orbit as expected, 7= -43±3 kms -1 , A"i=110±4 
kms -1 and A2=251±7 kms -1 . Using these values we estimate the pro- 
jected orbital semi-major axis and mass ratio as: asini=6.028±0.127 R© and 
q = ^=0.438±0.017. 

3.2. Effective temperature of the primary star 

The spectral type of a star can be derived from the multi-color photometry 
and spectroscopy, or both. Unfortunately, UBV measurements of the system 
are not available. Therefore we observed the system with the UBV stan- 
dard stars some nights and obtained the standard apparent visual magnitude 
and colors at mid-secondary eclipse as: V=9.231±5, U-B=0.101±15, and B- 
V=0.156±7 mag. The Q-paramet er of the Johnson's wide-b and photometry, 
defined as Q=(U-B) -0.72 (B-V) (| Johnson fc Morgan! fl95j ) . is independent 
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of interstellar reddening for early type stars. Using the observed (U-B) and 
(B-V) colors at secondary eclipse we calculated the Q-paramctcr as -0.011 
for the hotter star. Then, we esti mated the intrin s ic col or of (B-V)o = - 
0.021 mag from the tables given by iHovhannessianl (|2004l ). Since the light 
contribution of the secondary component does not exceed a few hundreths of 
magnitudes in the B and V passbands its contribution to the B-V color at 
maximum light may be ignored. This intrinsic color corresponds to an A0±1 
main-seque nce star. The average effective temperature deduced from the cal- 
i bratio ns o flBohm-Vitensel (119811) .Ide Jager fc Nieuwenhuiizenl (|1987l ). iFlower 
(jl996h and lDrilling fc Landolt] (fcoooh is about 9 810 K. The spectral-type un- 



certainty leads to an uncertainty of about 150 K in the effective temperature 
of the primary star. The difference between the observed and intrinsic B-V 
color, i.e. the interstellar red dening, is abou t 0.177 mag. 

KR Cyg was observed by iHilditchfc Hilll (|1975[ ) in the intermediate-band 
photometric system. One of their observations falls into the primary eclipse. 
Excluding this observation we calculated the colors of u — fo=1.555±0.011, 
b - t/=0.111±0.010 and c=0.970±0.022 mag. Taking E(b - y)=0.73E(B-V) we 
computed the intrinsic color of (b — y) =-0.018±0.010 mag which corresponds 
to an AO star. The infra-red magnitudes J=9.278± 0.022 H=9 169±0.022 and 
K=9.098±0.018 arc given in the 2MASS catalog JCutrill2003h . The infrared 
colors are calculated as J-H=0.109±0.031 and H-K=0.071±0.028 which cor- 
respond to a reddened B9±2 star. The spectral type of the primary star 
obtained both intermediate and infrared photometric systems is consistent 
what we derived from the UBV photometry 



3.3. Analyses of the light curves 

As we stated in §2.1 we have multi-color light curves of KRCyg obtained 
with the same instrumentation. We used the recent version of the ec lipsing 
binary light curve modelling algorithm of IWilson &: Devinnevl (|197ll) (with 



updates), as implemented in the PHOEBE code of Prsa & Zwitter (2005). 
The code needs some input parameters, such as the effective temperature of 
the hotter star and the mass-ratio of the system. These two parameters are 
the key-parameters for the solution which should be estimated before begin- 
ning the light curve modelling. The effective temperature of the hotter star, 
eclipsed in the primary minimum, has been estimated using the photometric 
indices as 9 810±150 K. The mass ratio, q=M^I Mi, is very important pa- 
rameter in the light curve analysis, because the WD code is based on Roche 
geometry which is sensitive to this quantity. The mass ratio of 0.438±0.017 
determined from the radial velocities was kept as a fix value. A preliminary 
estimate for the effective temperature of the cooler component is made using 
the depths of the eclipses. 

The logarithmic limb-darkening coefficients, x\ and X2, and bolometric 
limb-darkening coefficients, Xboi and yboi, were determined from tables by 
Van Hammc (1993) for the primary and secondary components, respectively, 
taking into account the effective temperatures and the wavelengths of the 
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observations. Standard values of bolometric albedos ( Rucinski 19731) . and 
the gravity-darkening coefficients (jLucvl 119671 ) for radiative and convective 
envelopes were used. The rotational velocities of the components are assumed 
to be synchronous with the orbital one. The adjustable parameters in the light 
curves fitting were the orbital inclination i, the surface potentials fix and f^, 
the effective temperature of secondary T2, the luminosity of the primary star 
Li. 

Using a trial-and-error method, we obtained a preliminary set of param- 
eters which represent the observed UBVR light curves. The code was set 
in Mode 2, a detached configuration, with coupling between luminosity and 
effective temperature and no constraints on the potentials. We used a sim- 
ple reflection treatment as MREF = 1 and NREF = 1 at the beginning of 
iterations. 

The iterations were carried out automatically until convergence and a solu- 
tion was defined as the set of parameters for which the differential corrections 
were smaller than the probable errors. After a few iterations convergence was 
obtained. The final results of the five-year multi-color light curves' simulta- 
neous analysis of KR Cyg are listed in Table 3. r\ and in this table are the 
mean fractional radii of the hotter and cooler components, respectively. The 
uncertainties assigned to the adjusted parameters are the internal errors pro- 
vided directly by the Wilson-Dcvinney code. We also ran models for KR Cyg 
in mode Mode-5, the secondary component filling its Roche lobe, in Mode-4, 
the primary filling its Roche lobe and in Mode-3, both components fill their 
Roche lobes. Convergence is obtained in Modes-4 and 3 but with larger sum 
of squared residuals. 

The computed light curves (continuous line) are compared with the obser- 
vations in four passbands in Figure 3. While there is a satisfactory fit between 
the observed and computed lights for the U passband the binary model does 
not represent the observed light curves successfully in the B, V, and R pass- 
bands. The residuals, in the sense observed minus calculated, are plotted in 
Fig. 4 for a better visibility. The differences between the model and the ob- 
served light curves are clearly seen at both maxima in three passbands and 
pa rticularly in the secondary eclipse in B and R passbands. As pointed out 
bv lShawi (|l990h all of the NCBs do not show O'Conncll effect, one maximum 
is higher than the other. If a NCB shows an asymmetry in the light curve, 
the first maximum is brighter than the second maximum. This difference is 
generally attributed to the secondary component, in which the visible surface 
at phase 0.25 is brighter than the other half. While the asymmetries reach to 
5% of the total light, the secondary components in these systems contribute 
less than 5% of the system's total light. Therefore, it is suggested that the 
visible hemispheres of the hotter components at the orbital phase of 0.75 are 
either cooler or obscured by a cool material. 

In the case of KR Cyg the light curve is almost symmetric in all passbands. 
At phases of about 0.15 and 0.85 the observed fluxes arc larger than those 
computed in the BVR passbands. In contrary, the observed fluxes are smaller 



KR CYGNI 



9 




Fig. 3. The observed five-year UBVR flux scale light curves and computed light 
curves for KRCyg. 
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TABLE 3 

FITTED PARAMETERS FOR THE FIVE-YEAR MULTI-COLOR LIGHT 

CURVES OF KRCYG. 



Parameter 


U 


B 


V 


R 


A (A) 


3500 


4300 


5500 


6700 


N ofas 


1752 


2825 


2804 


2825 


•'1 


0.620 


0.714 


0.611 


0.773 


X2 


0.871 


0.844 


0.773 


0.682 


9i 


1.0 


1.0 


1.0 


1.0 


92 


0.32 


0.32 


0.32 


0.32 


Ai 


1.0 


1.0 


1.0 


1.0 


A 2 


0.5 


0.5 


0.5 


0.5 


Ti (K) 


9810 


9810 


9810 


9810 


q=M 2 /Mx 


0.438 


0.438 


0.438 


0.438 


Li/(L 1+2 ) 


0.9791±0.0005 


0.9727±0.0004 


0.9467±0.0005 


0.9179±0.0006 



87.00±0.02 

T 2 (K) 5 580±10 

fii 2.8645±0.0012 
n 2 2.8111±0.0008 
n 0.4297±0.0002 
r 2 0.2989±0.0002 
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than those computed at the secondary eclipse in the B, V and R passbands. 
The components in the NCBs are nearly touch each other but their photo- 
sphere are not in physical contact. A gaseous stream and an energy exchange 
between the stars may be excluded. 

3.4. Search for the empirical derivation of the albedo 

The value of the local emergent bolometric flux and energy balance in 
a distorted star are significantly dependent upon the assumed values of th e 
gravity darkening exponent and the bolometric albedo faafert fc TwigJl98C)[) . 
The reflection effect is used to describe the mutual irradiation of the facing 
hemispheres of the components in a close binary system. For the stars of 
intermediate temperatures of about 10 000K the incoming radiation is almost 
completely absorbed by atoms and ions in the irradiated atmosphere which 
alters the temperature structure and causes an increased emergent flux. For 
the stars of convective envelopes half of the incoming r adiation is converted to 



bulk motions of the atmosphere gases (|Rucinskilll969f ). The theoretical value 



of bolometric albedo is assume d to be about 0.5 for th e stars having convec- 



tive atmospheres. Furthermore. IVaz fc Nordlundl (|1985l ) calculated bolometric 



reflection albedo for a particular reflecting main-sequence star with a temper- 
ature of 4 500 K and presented the results as a function o f angl e of incidence 



and relative incident flux. As pointed out by ICranmerl (|1993l ) the effect of 
mutual irradiation will be strongest on the line-of-centers points of the stars 
in a close binary system. The stars themselves iterate this effect between one 
another, e.g., hotter star heats up the cooler that in turn affects the illumi- 
nation back on the hotter star, and so on. The facing hemispheres of both 
components are heated up by multiple scattering process. 

Therefore we attempted to estimate the values of g-i and A2 from anal- 
yses of well-defined multi-color light curves of the system. WD code allows 
us to make parallel solutions by solving any number of subsets of adjustable 
parameters in the same computer run. By adopting the theoretical values of 
the less sensitive parameters of limb-darkening coefficients, gravity-darkening 
exponents and the limb-darkening coefficients we derived the basic set of pa- 
rameters i, Hi, ^2; L\, L2 and T2 from the analyses of multi-color light 
curves. The components of the system KRCyg are nearly contact and they 
rotate faster than 100 km s _1 , and have a large temperature difference of 
about 4 200 K. Therefore, the gravity-darkening exponent and the bolometric 
albedo for the components are expected to be different from the non-rotating 
and gravitationally undistorted stars. In the second step, we take gravity- 
darkening exponent, effective temperature, luminosity and effective albedo 
for the cooler component as adjustable parameters in the analysis of U, B, V, 
and R light curves separately. A few iterations indicate that there is no sig- 
nificant deviation in the gravity-darkening exponent between the theoretical 
and calculated values. The iterations are carried on excluding the gravity- 
darkening. After a few iterations a convergence is attained. The computed 
values of A2,and Ti are given with their standard deviations in Table 5. The 
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Fig. 4. Residuals, in light units, for the two solutions of KR Cyg's UBVR light curves 
from top to bottom: ^2=0. 5 (fixed, left panels) and A2 is an adjustable parameter 
(right panels). 



KR CYGNI 



13 



TABLE 4 

ABSOLUTE PARAMETERS OF KRCYG. 



Parameter 



Primary 



Secondary 



a (R ) 

Vy (km s- 1 ) 

q 

Mass (M ) 
Radius (Rq) 
log 9 ( C 9 S ) 

Teff (K) 

(vsin i) b s (km s _1 ) 
(vsin i) ca ic. (km s _1 ) 
log (L/L Q ) 
d (pc) 

U, B, V (mag) 
J, H, K s (mag) 



2.876±0.198 
2.594±0.055 
4.069±0.012 
9810±150 

141±7 

155±3 
1.750±0.032 

9.488±0.014, 
7.159±0.018, 



6.04±0.13 
-43±5 
0.438±0.017 

1.261±0.074 
1.804±0.038 
4.026±0.012 
5 580±150 
108±9 
108±2 
0.454±0.050 
411±12 
9.387±0.005, 9.231±0.005 
7.221±0.027, 7.249±0.027 



TABLE 5 

THE CALCULATED EFFECTIVE ALBEDOS AND EFFECTIVE 
TEMPERATURES FOR THE COOLER COMPONENTS OF KR CYC, 
AKCMI AND DO CAS. 



Star 


Parameter 


U 


B 


V 


R 


KRCyg 


A 2 


0.435±0.023 


0.225±0.016 


0.241±0.012 


0.304±0.012 


KRCyg 


T 2 


5600±27 


5985±14 


5653±12 


5336±12 


KRCyg 


L 2 


0.022 


0.041 


0.057 


0.070 


AK CMi 


A 2 


0.457±0.041 


0.132±0.038 


0.219±0.028 




AKCMi 


T 2 


5658±60 


5880±45 


5677±32 




DOCas 


A 2 


0.275±0.158 


0.377±0.127 


0.333±0.089 




DOCas 


T 2 


5441±147 


5615±105 


5427±74 





14 



SIPAHI ET AL. 



Cranmer 


^2 vaiuro 

(1993) and 


Claret ( 


J [JIGS 

2001) 



intrinsic and incident flux by 

<7T 4 = aT c 4 // + F„ 



■i dent 



(i) 



eff 



(2) 



where T e /y is the effective temperature of the star without external flux, 
Th is the effective temperature after irradiation, F r is the external radial 
flux and A is the bolometric albedo. The external flux F r is depended on 
apparent radius, effective temperature of the irradiating star, and the co- 
sine of the angle of incidence. The amount (l-A)F r is absorbed and the 
atmo sphere of the irra d iated star re-emit t h is exc ess energy later. More- 



over, iKirbivik &c Smith! (|1976l ) and iKirbivikl (|l982l ) proposed that not only 
the irradiated hemisphere of the star is heated in close binaries. The hori- 
zontal fluxes resulting from the incident radiation field can cause large cir- 
cular currents over the irradiated stellar surfaces and may penetrate into the 
non-directly irradiated hemisphere. As a result the irradiated star would in- 
crease its temperature and luminosity. The values given in Table 5 confirm 
this theoretical prediction. The empirically derived values of the A%, and 
T 2 from the light curves are plotted versus the wavelengths of the observa- 
tions in Fig. 5. The derived values of Ai are systematically smaller than their 
theoretical values of about 0.5. Moreover, the deviations appear to be de- 
pended upon the wavelengths. The largest deviation occurs about A4350. 
While the albedo gets smaller the effective temperature of the illuminated 
star increases. The computed value of albedo in the U-passband is close to 
the theoretically expected value for convective envelopes. The light contri- 
bution of the secondary star to the total light is too low for this passband. 
The influence of the mutual illumination in close binary systems 
is not only seen on the bolometric albedos but also in other prop- 
erties of the stars such as limb-darkening coefficients and gravity- 
darkening exponents. The reflection effect has a st rong influence on 
l imb-d a rkening coefficien t s, as a l ready noticed by 



Vaz &l Nordlund 



(|l985h . Nordlund &: Vaj (jl99dh . Iciaret fc GimenezJ (|l99oh . Later 
on, the effect of illumina tion on the limb-dark ening coefficients 
is studied numerically by Alencar &; Vazj (1999) using the Upp- 
sala Model Atmosphere code in convective line-blanketed atmo- 
spheres. Their results show that the limb-darkening coefficients 
of illuminated atmospheres are significantly different from the non- 
illuminated ones. They have concluded that the limb-darkening 
coefficients vary depending on the characteristics of external illu- 
mination such as the incidence angle, the amount of infalling flux, 
the te mperature of the illuminating star. In addition, 



Alencar et al. 



(1999) showed that external illumination changes the gravity-brightening 
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exponent of an illuminated star. They have proposed that the con- 
tact and semi-detached eclipsing binary systems should be better 
represented by illuminated atmospheres, since their components are 
closer to each other. Their numerical calculations show that the 
external illumination increases the values of gravity-brightening ex- 
ponent roughly in proportional to the amount of the incident flux. 
The classical value of 0.32 for the convective atmospheres may be 
too small for binary systems with close components. 

We also solved t he UBV light curve s of AK CMi, a similar system to 
KRCyg, obtained by Samec et al. ( 19981) . They estimated an effective tem- 
perature of 8 510K for the hotter star. We derived a value of 10 150 K using 
the UBV measurements of the system at four orbital phases, at maxima and 
in mid-eclipses. Applying the same procedure as in KRCyg we obtained A2, 
and T2 for the cooler component. They are listed in Table 5 and plotted in the 
middle panel of Fig. 5. Though the uncertainties are slightly larger both astro- 
physical parameters indicate very similar variations as in the cooler component 
of KRCyg. DO Cas is also a NCB and its UBV light curves were published 
by lOh fc Ahnl (|l992j ). We derived an effective temperature for the hotter 
star a s 8 700 K using the intermediate-band measurements by iHilditchfc Hill 
(119751 ) and the BVJHK mag nitudes given in the SIMBAD database. I nitial 
parameters were adopted from lSiwak. Zola &: Koziel-Wierzbowskal (|2010l ) who 
obtained spectroscopic mass-ratio and modelled the BV light curves. Their 
analysis resulted in near-contact configuration. Applying the same procedure 
we derive A2 and T2 for the irradiated star of DO Cas. The results are given 
in Table 5 and plotted in the bottom panels of Fig. 5. The uncertainties of 
the parameters are very large, originated from the large scatters in the ob- 
servational data, the effective albedos arc smaller than the expected value of 
0.5. 

In Fig. 6 we plot the observed U-B and B-V colors for the systems KRCyg, 
AK CMi and DO Cas. The bluer color especially in the U-B just in and around 
primary eclipse is clearly seen. About 50 per cent of the apparent disk of the 
hotter component of KR Cyg is obscured by the cooler star in mid-primary 
eclipse. We, therefore, suggest that multiple scattering of the light between 
the components of close binaries with large temperature differences results in 
bluer color, especially in the annular eclipse. 



4. CONCLUSIONS 

We obtained the five-year UBVR light curves and the radial velocities of 
the near-contact binary KRCyg. Analysing the multi-color data and mea- 
sured radial velocities we obtained, for the first time, the absolute parameters 
of the components. Comparison of the components' radii with the correspond- 
ing Roche lobes indicate that the system appears to be near-contact but not 
in contact. Combining the results of light-and radial- velocity analyses we 
obtained the absolute physical parameters of the system. Both components 
locate on the main-sequence of the Hertzsprung- Russell diagram. Systematic 
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Fig. 5. The empirical values of the effective albedos (left panels) and temperatures 
(right panels) are plotted against the wavelengths for the cooler components of KR 
Cyg (top panels), AK CMi (middle panels) and DO Cas (bottom panels). 




Fig. 6. The U-B and B-V colors against the orbital phase for KR Cyg (top panels), 
AK CMi (middle panels) and DO Cas (bottom panels). 
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behaviour of the residuals between the observed and fitted light curves is at- 
tributed to the effect of the mutual heating of the components, mainly in the 
case of the cooler star, of this near-contact system. Analyses indicate that 
the effective albedo and the effective temperature of the irradiated star are 
significantly altered. The empirically derived value of albedo is smaller than 
that expected for a convective star. While the albedo gets smaller the effective 
temperature increases. We also analyzed the UBV data, taken from literature, 
for the systems AK CMi and DO Cas. The changes of albedos versus the wave- 
length are very similar for AKCMi and KRCyg. The analysis of eclips- 
ing binary light curves involves a large number of parameters. The 
bolometric albedos, limb-darkening coefficients, gravity-brightening 
exponents are mostly kept fixed at the theoretical values. Never- 
theless, external illumination is significant in close eclipsing binary 
systems which heats the surface layers of the illuminated star, thus 
alters its physical properties. 
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